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Abstract 

We completely characterize the Narcowich-Wigner spectrum of Wigner functions associated with 
' pure states. 
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^ ■ 1 Introduction 
oo . 

, In the context of the Weyl- Wigner formulation of quantum mechanics, the derivation of criteria character- 
izing the phase space functions in terms of their quantum mechanical properties is a long standing problem. 
Such criteria should, at least, provide a simple and practical answer to the following questions: i) Which 
square integrable phase space functions are Wigner functions?; ii) Which Wigner functions are associated 
' with pure states?; iii) Which Wigner functions are associated with mixed states? and iv) Which Wigner 
functions are non-negative? 

Partial answers to these questions have been provided in the literature: Hudson [1] proved that Wigner 
functions of pure one dimensional states are everywhere non-negative iff the state is coherent. Soto and 
Claverie [2] generalized this result for higher dimensional systems. However, the analogous characterization 
of positive Wigner functions associated with mixed states remains an open question. Some results in this 
direction were obtained in [3] where the authors modeled the result of a measurement process on a system 
initially in a state described by a Wigner function F, in terms of the (non- negative) phase space distribution 
obtained by smoothing (convoluting) F with another Wigner function Fq. For fixed Fq^ this can be regarded 
as quantum dynamical map (a linear map that takes Wigner functions to Wigner functions) as long as the 
convolution F^-k F \s a, Wigner function for any Wigner function F. The authors then proved that this 
will happen provided Fq is point-wise non-negative. The procedure generates a large set of non-negative 
Wigner functions. On the other hand, in refs. [4, 5] necessary and sufficient conditions were derived for a 
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Gaussian phase space function to be a Wigner function. Gaussians are obviously important for their role 
in optics [6] and because they arc bona fide probability measures. Moreover, they provide kernels Fq for 
quantum dynamical maps in the above sense. Strictly speaking, Gracia-Bondfa and Varilly [4] generalized 
the result of [3] by proving that: (i) if Fq is any point-wise non-negative phase space function (Wigner 
function or not), then F^-k F is a Wigner function for any Wigner function F, and (ii) if Fq and F are 
Wigner functions, then Fq-mF is always non-negative (albeit possibly not a Wigner function). Using these 
results, they constructed explicit examples that preclude the Hudson, Soto, Claverie theorem for mixed 
states; i.e. there are non-Gaussian positive Wigner functions associated with mixed states. 

Other authors considered whether certain criteria (namely the uncertainty principle) would be useful in 
assessing whether a phase space function is a Wigner function. In [7, 8] it was shown that the uncertainty 
principle does not determine the quantum state. Conditions for Gaussians in phase space to be Wigner 
distributions were expressed in terms of the so-called symplectic capacity of the associated Wigner ellipsoid 
[9]. In [10], the multidimensional Hardy uncertainty principle was extended and expressed in the context 
of Wigner quasi-distributions. 

The concept of the Narcowich- Wigner (NW) spectrum was introduced and explored in refs.[ll], [12]. 
Its main purpose was to provide a criterion to classify certain phase space functions and a framework 
suitable to generate Wigner functions displaying specific properties (namely being non- negative) . With 
each square intcgrable phase space function F, wc may associate its NW spectrum, denoted as W(-F), 
which is a compact subset of IR. If this subset contains Planck's constant h, then the phase space function 
F is a Wigner function and represents some quantum mechanical state. If, on the other hand, G W(F), 
then F is a point-wise non-negative function, that is a classical probability measure. We see that, quite 
straightforwardly, the NW spectrum provides a criterion for questions i) and iv) above. 

On the other hand, the results of [3], [5] were all reformulated and unified in the context of NW spectra 
in ref.[12]. Various authors [11], [13] also speculated about the possibility of using the NW spectrum to 
generate the entire subset of Wigner functions which are everywhere non-negative (i.e. contain in their 
NW spectra). While this has still not been fully achieved, by resorting to the NW spectrum, Narcowich 
[11] proved that one can generate a large set of positive Wigner functions of mixed states by convoluting 
suitable Wigner functions. This was basically a generalization of the previous results of [3] and [4]. Other 
potentially interesting applications of the NW spectrum are in the field of signal processing [14] and in the 
study of the classical limit of quantum mechanics [15]. 

An unresolved problem is that of using the NW spectrum to characterize pure and mixed states and 
thus to answer the questions ii) and iii) above. In its general form this is still an open issue. For Gaussian 
states, however, the spectrum is of the form W{F) = [— r/o, Vo] (eq-(20)) and it is known that ryo = ^ iff the 
state is pure and ijq > h iS it is mixed (770 < fi iS F is not a Wigner function). 

In this work we will prove an important result towards the solution of the above problem: we will 
complete the characterization of the NW spectrum of pure states. In fact the purpose of this paper is to 
prove the following theorem. 

Main Theorem Let F G L^(]R^'', dS_) be the Wigner function of a pure state. 

(i) If F is a Gaussian, then W(F) = [—%,%]. 

(ii) If F is non-Gaussian, then W{F) = {—h,h}. 

Notice that, as mentioned before, the result (i) in the Main Theorem is already known [11]. We include it 
here for completeness. 

We will also prove an auxilary result concerning the zeros of the so-called Husimi function [16]. This 
result is stated in proposition 2.4. 

The Main Theorem is of course consistent with the Hudson, Soto, Claverie theorem: Gaussian states are 
the only pure states that contain zero in their NW spectrum. Unfortunately, the analogous characterization 
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of the NW spectrum of mixed states remains an open issue. As shown in [13], mixed states can have a 
very intricate Wigncr spectrum: intervals, sequences, or combinations of both. 

This paper is organized as follows: Section 2 introduces the notation, the main definitions and reviews 
the fundamental properties of the NW spectrum. In section 3 we prove our main results concerning the 
NW spectrum of pure states. We start by proving some auxiliary results for the Husimi function and review 
some properties of Bargmann transforms in section 3.1, and in section 3.2 we prove the Main Theorem. 
Finally, in section 4 we discuss some generalizations and possible applications of these results. 



2 The Narcowich-Wigner spectrum 

Let us first settle the preliminaries. Wc shall consider a d-dimensional system on a flat phase space 
T*M ~ (IR'^)* X IR"^ ~ IR^'^ with a global Darboux chart ^ = {q,p) G K^'' in terms of which the symplectic 
form reads: 

Y -i-dxd ^dxd J 

Here the superscript "T" denotes matrix transposition. We shall use the standard norms in real, ||a;||]R" = 



xf + ■ ■ ■ + x"^, and in complex vector spaces ||^;||([;" = \/\zi\'^ + • • • + \zn\'^- To avoid a proliferation of 
subscripts, we shall just write or as it will be clear from the context the vector space we are 
referring to as well as its dimensionality. As usual L^{JR",dx) denotes the space of square integrable 
complex- valued functions on IR" with respect to the Lebesgue measure. We may define in L^{JR",dx) the 
inner product: 

<il;\(f)>= [ ^^(j){x)dx, ^p,(f)e L'^{WC,dx). (2) 

We shall also consider the vector space T of Hilbert-Schmidt operators acting on the Hilbert space Ti. = 
L^('Mf',dq) of our system, and which admit a kernel representation of the form: 

m ^ {A^){q)=J^,A{q,q')^{q')dq', 
with A(q,q') G L'^{lR?'^,dqdq'). The Weyl-Wigner transform is the isomorphism [17]: 
Wr,: T ^ L\M^^,dC) 

where rj is some real, positive constant. The Weyl-Wigner transform may be extended to other spaces, but 
this is all that will be necessary in this work. For a state il^{q) G H, we define the corresponding density 
matrix to be the operator G T: 

{p^<P){q)= / r^{q,q')(f){q')dq', (5) 



with r^{q,q') = i^{q)ip{q')- The 77-Wigner function associated with V' is [18]: 

wMm) = j^d^vipM) = e-'^^-y/^W^m + y)dy. (6) 

If r] is equal to Planck's constant we shall denote Wfi{ijj, simply by Wigner function or Wigner 

quasi-distribution in accordance with the literature. Moreover, it is straightforward to prove the following 
useful identity, also known as the Moyal identity [19]: 

J^^^Wr,{^,m)w,{ci^,m)d^ = < ^\4> > (7) 
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where Wr^{ilj,ilj),Wr]{(f),(f)) are ry-Wigner functions associated with states il'jcf) E L^{lR'^,dq). 
We may consider statistical mixtures of pure states: 

P = ^VaP^ac^ Pa > 0, = 1, (8) 

a a 

which yield, via the Weyl-Wigner transform, the mixed state ry-Wigner functions: 

Wr,{pm = E^'«^'?(^a,^a)(6- (9) 
a 

The Wigner functions thus represent the states in the quantum phase space. An important issue is the 
identification of necessary and sufficient conditions for a phase space function to be a ry- Wigner function 
[20]. A well-known necessary condition is stated in the following proposition [20]. 

Proposition 1.1 If is a r]-Wigner function, then: 

In addition, the equality holds if and only if F'^ is associated with a pure state. 

The integral on the left-hand side of the previous inequality is called the purity of the system. 

The concept of NW spectrum was introduced in ref.[ll]. It appears naturally if one wants to formulate 
necessary and sufficient conditions (the KLM conditions) for a phase space function to be a Wigner function 
or, more generally, a 77- Wigner function. In addition, we are able to (partially) characterize different types of 
Wigner functions in terms of their NW spectra, namely the non-negative Wigner functions, the Gaussians, 
and pure or mixed states [12]. For this purpose it is useful to define the symplectic Fourier transform. 
Let G(^) G L^(IR^'^, c/^). We define its symplectic Fourier transform according to: 

: L2(1R2°') ^ L2(IR2^) : G(0 ^ G(a) = (.FG)(a) = / G(0 exp (iO(e, a)) (11) 

The inverse transform is given by the formula: 

G{0 = (^-^G)(0 = j^^^ G{a) exp a)) da. (12) 

Definition 1.2 The symplectic Fourier transform F{a) is said to be of the a-positive type (a G IR) if 

the m X m matrix with entries 

Mjk = F{aj - ajk)exp [^0(0^,0^)) (13) 

is hermitian and non-negative for any positive integer m and any set of m points oi, • • • , Om in the dual 
of the phase space. By abuse of language, we shall often say that a phase space function -F(4) is of the 
a-positive type, by which we mean that its symplectic Fourier transform F{a) is of the a-positive type. The 
Narcowich- Wigner spectrum of a phase space function F{^) G L^(]R^'^, c/^) is the set: 

F{a) is of the a-positive typej . (14) 

With this definition we may now state the KLM (Kastler, Loupias, Miracle-Sole [21], [22], [23]) conditions. 
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(10) 



W{F) = {a G IR 



Theorem 1.3 The phase space function F{$,) is a rj-Wigner function (rj > 0), iff its symplectic Fourier 
transform F{a) satisfies the KLM conditions: 

(i) m =1 (15) 

(a) F{a) is continuous and r) G VV(-F). (16) 



This theorem is a twisted generahzation of Bochner's theorem, which states that non-negative functions 
(i.e. classical probability densities) are those of 0-positive type. To proceed we introduce the convolution 
ofF,Ge L2(IR2d^ ^hich is defined by: 

{F*Gm= F{C-^')G{C')dC', (17) 

and is continuous, bounded and vanishes for ||^|| oo. Moreover, the symplectic Fourier transform of the 
convolution is tantamount to point-wise multiplication: 

iT{F * G))(a) = (Hma) ■ {HG))ia). (18) 

Let us now recapitulate some of the properties of the NW spectrum. 

Proposition 1.4 Let F,G e L'^ {JR'^'^ , d^ ■ 

(i) If f is continuous, then W(F) is hounded and closed 

(ii) a G W{F) ^ {-a) G W{F) 

(Hi) {a + a'l a G >V(F), a' G W{G)} C W{F * G), where {F * G) is the convolution in (17). 

The proof of these properties can be found in [13]. Counter-examples to the reciprocal of property (iii) are 
also given in rcf.[13]. In view of property (ii) of proposition 1.4, we may assume without loss of generality 
that a > 0. The results of [4], [5] for Gaussian states can be expressed in terms of the NW spectrum as 
follows (see [11] for a proof): 

Lemma 1.5 Let F be a Gaussian 



= ^^15^ exp [-(^ - ^o)^A(C - Co)] , (19) 

where A is a real, symmetric, positive definite 2d x 2d matrix and ^0 € M^'^. Then the NW spectrum of F 
is of the form: 

W(F) = [-770,770], (20) 
where rjo = max(W{F)). That this maximum exists is guaranteed by property (i) of proposition 1.4- 

The following lemma states the well known results of [3], [4] mentioned in the introduction. We include it 
here to shown the usefulness of the NW spectra. Indeed, using the previous results, the proof is trivial, as 
shown in [12]. 

Lemma 1.6 Let F, Fq be Wigner functions. 

(i) The convolution Fq * F is point-wise non-negative. 

(ii) If Fq is such that G W{Fq) or 2h G yV{Fo), then the convolution Fq* F is a Wigner function. 
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Proof If F,Fo are Wigner functions, then {—h,h} C W{Fq) D W{F) and thus, from proposition 1.4 (iii) 
G W(-Fo * F), which means that FQ-k F is point-wise non-negative. On the other hand, if Fq is such that 
{0,±n} C W{Fo) or {±n,±2n} C >V(Fo), then from proposition 1.4 (iii), {-^,0,^} C >V(Fo * F), and 
thus Fo-k F is a Wigner function and point-wise non-negative. □ 

This is a powerful way of generating positive Wigner functions, notwithstanding the fact that the method 
does not exhaust the entire set of positive Wigner functions. 

3 Main results 

3.1 Husimi function and Bargmann transform 

We start by deriving some resuhs which wiU be useful to prove the Main Theorem. The following restricted 
version of Hadamard's theorem for functions of several complex variables was proved in [2]: 

Theorem 2.1 (Hadamard, Soto, Claverie) If F{z) is an entire function on C", with order of growth 
p and without zeroes, we have: 

F{z) = e^p{P{z)) (21) 

where P{z) is a polynomial of degree r < p. 

We recall that the order of growth p of F{z) is given by: 

log(logM(i.)) sup |F(.)| 

R^oo logi? ||2;||=iZ ^ ' 

Let us now consider the coherent state 

Mq) = (^^)~^ exp (^-M! + ^ . ^ _ ^\\Rez\f^ , zeC''. (23) 

Moreover, let W^ii/Jz, ^2) and Wn{il^, i/j) denote the Wigner functions associated with ipz and ip G L^{JR'^, dq), 
respectively. 

Definition 2.2 Let ^ G L^{lR'^,dq). The convolution of the form 

QiO = (Wnii^z, i^z) * Wnii^, V))(0 (24) 
is called a Husimi function. The Bargmann transform of is defined by [17], [24]: 

F(z)=<V|V. >e^"''''"' = -^_/^,exp||-M! + ^.gj^^g, ^eC'^. (25) 

The following lemma provides an important characterization of the Bargmann transform: 

Lemma 2.3 Let tp G L^{M!^, dq) and let F{z) he its Bargmann transform. Then F(z) is an entire function 
on C^. Moreover, the order of growth p of F{z) is at most 2. 

Proof Since the integral in (25) converges uniformly in any compact subset of C^, we conclude that 
F{z) is an entire function on C^. On the other hand, since ip^ipz are both normalized, we have from the 
Cauchy-Schwarz inequality: 

\F{z)\ = elll«-ll'| < V'lV. > I < etil^^^ll'. (26) 
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Consequently, M{R) < e^^ and thus p < 2. a 

The following result will be useful to prove the Main Theorem. 

Proposition 2.4 The Husimi function Q{^) has no zeroes iff the state ip G L^(JR'^,dq) is a Gaussian. 
Proof If we compute Wji{'tpz,'4^z) explicitly (cf.(6,23)), we obtain: 

Wn{A,i^zm = (^e^P (-^11^ - ClP) (27) 

where 

( = [x, k) = {h Re{z),n Im{z)). (28) 

The convolution (24) thus reads: 

We conclude that if is a Gaussian, then Wn{il^,il^) is a Gaussian and so is Q. Then, evidently, Q has no 
zeroes. 

Conversely, let us assume that is not a Gaussian, but that it has no zeroes. This means that the 
Bargmann transform F{z) also has no zeroes. Indeed, we have: 

\F{z)\^ = e'^ll^-ll'l < ^IV-. > P = {27rnfe''\\^'^'\\' J^^a Wn{tPz, ^z^OWniA m)dC = 

(30) 

= 2<^e'^ll^-ll' 1^2, exp (-i||C - e'lP) Wnii^,m')de = (27r?i)'^e^ll«-ll'Q(2C), 

where we used (7). This means that Q has zeroes iff F(z) has zeroes. From the Hadamard, Soto Claverie 
theorem (Theorem 2.1), since F{z) is an entire function on C^, with order of growth p < 2 and, by 
hypothesis, without zeroes, then it must be a Gaussian. We then have, in particular, that: 

F{W) = -^ f ,^w(--^ + iyq]i^dq, y G IR^ (31) 
(7r^)4 Jr'* \ 2n J 



is a Gaussian. As F{iy) is the Fourier transform of e tlj{q), this is only possible if e tjj[q) is also 
a Gaussian. Hence ilj{q) is a Gaussian, which contradicts our assumption.n 

3.2 Proof of the Main Theorem 

The following lemma restricts the NW spectra of pure states. 
Lemma 2.5 Let F be the Wigner function of a pure state. Then: 

W{F) C [-n, n] (32) 

Proof From proposition 1.1, we know that, since _F is a Wigner function the equality in (10) holds for 
r] = h. Moreover, if r/ > is some other element in W(-F), then again from (10), we know that: 
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It follows that < rj < h. The lemma is then an immediate consequence of property (ii) of proposition 
1A.O 

The main theorem of this work refines the previous lemma. 
Proof of the Main Theorem 

(i) This is well-known. Wc prove it here for completeness. Let F be a Gaussian (19). From lemma 1.5, we 
know that for F to be a Wigner function (pure or mixed), then the inclusion must hold: 

[-n,n]c[-no,vo]=w{F). (34) 

But from lemma 2.5, if F is associated with a pure state, then the inverse inclusion must also hold. We 
conclude that: 

w{F) = [-n, n] (35) 

(ii) Let us now assume that F{^) is the Wigner function of a pure state, but that it is not a Gaussian. 
We already know from lemma 2.5, that its NW spectrum is contained in [—h, h]. Let us then assume that 
there exists ry G ] 0, ^ [ PI W{F). Let us consider the family of Gaussians: 

GaiO = -ATae^p(--\m'), «>0 (36) 

The NW spectrum is W(Gq) = [—a, a] [11]. The convolution of two such Gaussians is again a Gaussian 
in this family. In fact, the convolution verifies the semi-group law: 

(G„ * Gfi\i) = Ga+piO, > (37) 

Let us then define the function: 

QiO = (Gn-^ * Gr, * Fm = {Gn * Fm = (Gn-r, * QoW, (38) 

where (5o(0 = (Gr^ * f)iO- The Wigner spectrum of Grj is yV(G^) = [— ??,??], whereas that of F contains 
the set {—h, —Tj, Tj, Ti} by assumption. Prom property (iii) of proposition 1.4, we conclude that {—Ti, 0, 1%} C 
W(Qo)- This means that Qq is both a Wigner function and, according to Bochner's theorem, everywhere 
non-negative. Since Q{C) = {Gn*F){^) is the convolution of a Gaussian of width fi with a Wigner function, 
it is a Husimi function. Moreover, since F is not a Gaussian, then Q{C) least one zero (proposition 

2.4): 

Qiio) = (39) 



Prom (38), we then have: 



This is only possible, if Qo(0 ^ negative part. However, that is contradictory with the fact that Qq is 
of 0-positive type. This proves that there can be no 77 € ] 0, /i [ H W{F). Altogether, this means that the 
NW spectrum of a non-Gaussian pure state is contained in {—h, 0, h}. However, since non-Gaussian pure 
states cannot have non-negative Wigner functions (by the Hudson, Soto, Claverie Theorem), we conclude 
that ^ W{F) and thus W{F) = {-h, hj.n 
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4 Concluding Remarks 



A few remarks are now in order: 

• The Main Theorem admits a trivial generaHzation. We can easily replace h by any 77 > 0. And thus, if ijj 
is a Gaussian, then the 77-Wigner function W^('0,'0)(O Wigner spectrum [—rj,rj], otherwise its Wigner 
spectrum is {—r],r]}. 

• The converse result of the Main Theorem remains as an open question. That is: suppose F G L^(IR^'^, dS,) 
is such that F{a) is continuous and has NW spectrum W{F) = {—h,h}. Does that mean that F is the 
Wigner function associated with a non-Gaussian pure state? Or could it be a mixed state? If we could 
prove this result then pure and mixed states could be distinguished exclusively in terms of their NW 
spectra. 

• In [11] Narcowich stated the following conjecture: 

The convolution of a fixed Wigner function F^ with any other Wigner function Fj^ is again a Wigner 
function iff the NW spectrum of F^ contains or ±2^. 

As we mentioned before, the sufficiency of the conditions is trivial to prove (see lemma 1.6). Narcowich ar- 
gued that the necessity is proved, if we assume that (i) W{F*G) = W{F)+W{G) = {a + a'\a € W{F), a' G >V(G)} 
for any F,G€ L^(]R^'^, d^); and (ii) there exists a family of Wigner functions Ff^(^) defined for all suf- 
ficiently small e > and satisfying W{F^) C [—h — e, —h + e] U [h — e,h + e]. Let us reproduce here his 
argument for completeness. If F^ is a Wigner function and F° * Fji is yet another Wigner function for any 
Wigner function Fn, then h G W{F^*Fn). From hypothesis (i), it follows that h G W(F^) + >V(F|). If we 
take e \ 0, wc get h G yV(F^) + { — h, %}, which is possible if and only if cither or ±2% belong to W{F^J. 
The problem with this argument is that, as mentioned before, hypothesis (i) is not true in general [13]. 
However, that does not mean that the conjectured theorem is wrong. To the best of our knowledge no 
counter-example has ever been found. The problem might be that the hypotheses suggested may be too 
restrictive. Let us just assume that for the fixed Wigner function F^, we can always find a non-Gaussian 
pure state V such that W(F^ * WnH'ji')) = W{F^) + W {WnH^ , ■ Notice that we do not require this 
to be true for all Wigner functions. We only need this to be true for one state Wti{ip,il)). We then have, 
from the Main Theorem: h G W(F^ * Wfi{tlj,4')) = W(-F^) + {-h, +h}. And the necessity would thus be 
proved. The main result of this paper shows that at least hypothesis (ii) is easily met. 

• Our lemma 2.5 has an important consequence in the context of open (dissipative quantum systems). In 
this framework one usually regards the purity (10) of the system interacting with an external environment 
[25] as an indication of dissipative effects. In fact it tends to decrease. In most cases if -Fk(0 is the Wigner 
function of the system at the initial time, then the Wigner function will usually take the form Gt{^) * Fji{^) 
at a later time t (up to composition with a linear symplectic transformation, which leaves the NW spectrum 
unchanged [11]). Here Gt(^) is the Green function of the master equation. If the Narcowich conjecture 
mentioned above is correct, then Gt should be such that ±2h G W{Gt) or, alternatively that G W(Gt). 
Moreover, if Fji is associated with a pure state (F^ = Wn{il^, then from lemma 2.5, we infer that Gt*F% 
must be associated with a mixed state, whenever there exists 77 7^ in W(Gt). 
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